Exact dynamics of quantum correlations of two qubits coupled to bosonic baths 
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The counter rotating-wave term (CRT) effects from the system-bath coherence on the dynamics of 
quantum correlation of two qubits in two independent baths and a common bath are systematically 
investigated. The hierarchy approach is extended to solve the relevant spin boson models with 
the Lorentzian spectrum, the exact dynamics for the quantum entanglement and quantum discord 
(QD) are evaluated, and the comparisons with previous ones under the rotating-wave approximation 
are performed. For the two independent baths, beyond the weak system-bath coupling, the CRT 
essentially changes the evolution of both entanglement and QD. With the increase of the coupling, 
the revival of the entanglement is suppressed dramatically and finally disappears, and the QD 
becomes smaller monotonically. For the common bath, the entanglement is also suppressed by the 
CRT, but the QD shows quite different behaviors, if initiated from the correlated Bell states. In the 
non-Markovian regime, the QD is almost not influenced by the CRT and generally finite in the long 
time evolution at any coupling, while in the Markovian regime, it is significantly enhanced with the 
strong coupling. 
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I. INTRODUCTION 

Quantum correlation, originated from the coherent su- 
perposition in quantum mechanics P, [H, has been at- 
tracting persistent attention ranging from quantum infor- 
mation, quantum optics, to many-body physics 0-111 . As 
the characteristic trait of the quantum correlation pro- 
posed by Schrodinger @ without the classical counter- 
part, quantum entanglement played the crucial role in 
the highly efficient quantum computation and quantum 
information processing 0, Q ■ I* has also been applied to 
identify the phase transition in quantum many-particle 
systems Q. However, there exist exceptions where the 
entanglement seems unnecessary, such as in the Grover 
search algorithm (lol . [ll[ and the deterministic quantum 
computation with one pure qubit (DQCl)[l3|- It is evi- 
dent that other kind of the nonclassical quantum corre- 
lation rather entanglement dominates the computational 
implementation of DQCl 13]. Zurek and Vedral de- 
scribed one particular measure to quantify all nonclassi- 
cal correlation, termed as quantum discord (QD) il^ - [l6| . 
Only when the QD disappears, the information can be 
safely obtained by locally measuring the distantly sepa- 
rate subsystem without disturbing the bipartite quantum 
system, where the corresponding state is fully classical. 
Recently, the QD has successfully been used to study the 
quantum phase transition of the critical system s I17h22| , 
operational meanings for quantum processors [23|. and 
state preparation [24| . 

It is known that the realistic system inevitably in- 
teracts with the environment, resulting in the decoher- 
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ence, which is of fundamental importance in the quan- 
tum information processin g an d measurement [25l [26| . 
Specially, Yu and Eberly [271 found that the dynami- 
cal behavior of the entanglement is significantly differ- 
ent from the single qubit in a pair of qubits separately 
coupling with Markovian baths, which can be modeled 
from the well-known spin boson models [28]. It shows 
finite time disent ang lement called as entanglement sud- 
den death (ESD) [23] ■ However, this was investigated un- 
der quantum master equation with the Born and Markov 
approximation, which assumes that system-bath interac- 
tion is weak and the relaxation of the bath is much faster 
than that of the system. Then further works have been 
done in the non-Markovian regime in the framework of 
the rotating-wave approximation (RWA) to reveal novel 
characters of the quantum correlations compared to the 
Markovian effect }29n34] . which is only valid under the 
condition of weak system-bath interaction. Besides, en- 
tanglement dynamics beyond the weak coupling spin bo- 
son models has also been exploited based on the opti- 
mal polaron transformation [35| , where the higher order 
terms are neglected. It is pointed out that the corre- 
sponding accuracy in the wide parameter zone should be 
carefully analyzed [s^. Therefore, It is highly desirable 
to solve the dynamics without performing any approxi- 
mation. 

Beyond both the Born-Markov approximation and 
RWA, Tanimura et al. Q have developed an efficient 
hierarchy approach, which has been later extensively ap- 
plied to some chemical and biophysical systems [38, d9] . 
Recently, this method has also been extended to study 
the dynamics of entanglement for two qubits coupled to 
a common bath [4^ |41[. However, the dynamics of QD, 
one recently mostly studied quantum correlations has not 
been investigated beyond the Born-Markov approxima- 
tion and the RWA. Moreover, the exact study of the dy- 
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namics of both entanglement and QD for two qubits cou- 
pled to their own baths are not available in the literature, 
to the best of our knowledge. This case is designed for 
two remote qubits each interacting locally with its own 
environment. Such a physical condition is also relevant 
to the quantum information science and quantum com- 
putation. 

In the present paper, to give a comprehensive picture 
of the two qubits coupled to their own bath and com- 
mon bath, we will extend the hierarchy equation to these 
two kinds of baths to study the dynamics of the pairwise 
entanglement and the QD in both Markovian and non- 
Markovian regimes. The effect of the counter-rotating 
term on these dynamics in a wide coupling regime will 
be systematically studied. This paper is organized as fol- 
lows. In Sec. II we briefly review the definition of the 
QD. In Sec. Ill we exactly solve the reduced qubits den- 
sity matrix of two independent spin boson models, and 
investigate its quantum correlation compared with that 
in RWA. In Sec. IV the quantum correlation of two sepa- 
rate qubits interacting with a common bath is analyzed. 
Finally, a summary is given in Sec. V. 



II. QUANTUM DISCORD 

The QD is one main route to fully measure the non- 
classical correlation, which can be extended from the 
classical information theory 0, [T^ - [l6| . It is interpreted 
as the difference of two quantum mutual correlations of 
subsystems A and B, before and after local measure- 
ment operated on subsystem B [3 ■ The total quan- 
tum correlation of two subsystems is determined by their 
joint density matrix p as X{p) = S{pa) + S{pb) — S{p), 
where Pj\(b) = ^t^b{A){p}- The von Neumann entropy 
S{pa) — — Tr{pa Inpa}. The other quantum mutual in- 
formation J'(p), which is equivalent to in classical 
information frame, is derived by locally measuring sub- 
system B with a complete set of orthonormal projectors 
{n^}, where k is one outcome state of B. After this 
measurement, the joint density is reduced to conditional 
counterpart as pk = [(I®nf )p(I®nf where the 

corresponding probability is pk = Tr{(I(X)n^)p}. Then 
the mutual information based on specific performance is 
shown as Q = S{pa) — J2kPkSipk)- Hence, we should 
maximize it to capture all classical correlation by obtain- 
ing Jip) = max|nB}{Q}- Finally, the QD is defined as 



V{p)=I{p)~Jip). 



(1) 



The quantum nature of the system can be explicitly ob- 
served from the QD. It is only zero in pure classical cor- 
relation condition, and has the same values as the en- 
tanglement in pure state. While it remains finite even in 
separable mixed states, where the corresponding entan- 
glement may completely disappear [14 |. 



III. TWO INDEPENDENT SPIN-BOSON 
MODELS 

The Hamiltonian of two independent qubits coupled 
two independent bosonic baths A and B is given by if = 
Hs + Hb + Hsb, reads 



=A.B 



a=A,B;k 

Hsb = '^^.(gtKk + gTbU). (2) 

a=A,B;k 

where a — A, B , ct^ (/3 — x, y, z) is the Pauli operator 
with Zeeman energy uja of spin /3, h\ ^. {ba,k) creates (an- 
nihilates) one boson with frequency in bath a, and 
is the coupling strength between the system and the 
bath. 

The spectra density of the bosonic bath is selected as 
the following Lorentz distribution 



(3) 



which can describe an bosonic field inside an imperfect 
cavity mode with the system-bath coupling strength 
Xa- The corresponding correlation function at zero tem- 
perature is shown as (jA13p 



Caa(w) = exp[-(7a + 



(4) 



Then by applying the Feynman- Vernon influence func- 
tional procedure in Appendix, the hierarchy equation is 
derived by (IA16I) 



d_ 

dt 



(5) 



-[iHg + n- pA+rn- PB)Pn,rni^) 



k=l,2 

Xaia 



Y nkm + i-lf+'Q'^]pn^e,Mt) 



k=l,2 

+ Yll {'~^)^Q^Pn,'m+e, 
k=l,2 

, XbIb 



Y ^^\Qb + {-lf+^Ql]Pn.n.-eS)- 



fc=l,2 

The commutation relations B — AB—BA and A°B = 
AB + BA. The unit vectors ei = (1,0), 62 = (0,1), 
and fi — (?ii,n2), rn — (7711,7712). Qa (a = ^1 B) in 
the system-bath interaction represents cr°. At the initial 
state, only P(o,o).(o,o)(0) is not zero. For the simplest con- 
dition in the following calculations, the two subsystems 
are identical as = wo, Aa = A, 7a = 7. 
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For comparison, we thereby briefly review the main 
results derived in the RWA, where system-bath interac- 

tion becomes H^, = EaAdk^^aM + gT^-blk)- When 
the initial system reduced density matrix in the basis 
{|1) = 111), |2) = |10), |3) = |01), |4) = |00)} is X form 



Pab(O) = 



then we have 



^Pll(O) pi4(0)N 

P22(0) P23(0) 

P32(0) P33(0) 

^P4l(0) P44(0)y 



(6) 



(a) 

^- ?.=0.02,t=5 Hierarchy 

I 'l ?.=0.02,t=5 RWA 

1 i ?.=0.02,f=0.1 Hierarchy 

I \ ?.=0.02,t=0.1 RWA 

I I 


\ (b) 

1 •- ?.=0.5,f=5 Hierarchy 

' \-. ?.=0.5,f=5 RWA 

1 V ?.=0.5,f=0.1 Hierarchy 

> V. ?.=0.5,f=0.1 RWA 

1 V 


, (c) 

- ?.=2,f=5 Hierarchy 

1 :- X=2,f=5 RWA 

1 ?.=2,f=0.1 Hierarchy 

' i ' ?.=2,f=0.1 RWA 

1 1 

1 J . 


\. \ ../ 


' V- 

\ \ 


1 » 

\ \ 1 
^ x ' 



Pll(t) - PlMPt. P22{t)= P22{^)Pt+ PlMPt{l-Pt). |«-^ 

P33(i) = P33(0)Pt+Pll(0)Pt(l-Pt), Oo.2 

P44(i) = 1 - [Pll(i) - P22(i) - P33(<)]> 

Pl4(i) = Pl4(0)P(t), P23W =P23(0)P(t), (7) 
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(d) 




(e) 




(f) 


?.=0.02,f=5 Hierarchy 




?.=0.5,f=5 Hierarchy 




?L=2,f=5 Hierarchy 


?.=0.02,t=5 RWA 




?.=0.5,f=5 RWA 




?.=2,f=5 RWA 


■\ X=0.02,t=0.1 Hierarchy 




J.=0.5,f=0.1 Hierarchy 




A.=2,f=0.1 Hierarchy 
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1 V 
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where 

Pt = e-'<'[cos{Rt) + sin(i?0]' (8) 



with R — A7/2 — 7^/4. It is known that the bath 
relaxation time is T;,sal/7, and the system-bath correla- 
tion time is rr~l/A. From the time dependence of the 
function Pj, we immediately know that the Markovian 
regime corresponds to 7/2 > A where Pi decays expo- 
nentially, while the non-Markovian regime to 7/2 < A 
where Pt shows oscillating behavior describing the coher- 
ent process between the system and the bath. This non- 
Markovian condition has been realized in cavity quantum 
electrodynamics having Rydbcrg atoms inside the Fabry- 
Perot cavities with 7/A«0.1 |42l- We should point out 
here that the Markovian and non-Markovian regimes are 
only defined in the framework of the RWA, which can 
not be well defined in the exact treatment without the 
RWA. Throughout the present paper, we still use such 
notations. Doing so is only for comparisons and analysis. 

As usual, the initial Bell states with anti-correlated 
spins and correlated spins will be considered in the 
present paper 

|$(a)) = a|10> + Vl -a^lOl), (9) 
\^{a)) = a|00) + -a2|ll>, (10) 

which both obey X form in Eq. (|6]) during the evolu- 
tion. So it is straightforward to derive the concurrence, 
a pairwise entanglement [29j 

C^(t) = niax{0,2av/l -a^pj^ (H) 

C*(i) = niax{0, 2VT~a2 Pt[a - VT~~a2(l - Ft)]}- 

(12) 

The analytical expression for QD in this case is quite 
complicated and also only consists of a and Pf, which is 
implied in Refs. [H IH 



FIG. 1: Dynamics of quantum correlation of two qubits cou- 
pled to two independent baths with initial anti-correlated Bell 
state i$(-i=)>. a;o = l, 7 = /A- 



It is shown analytically above that the dynamics of 
both entanglement and QD under the RWA can be dis- 
tinguished in the non-Markovian and Markovian regimes, 
characterized by the ratio 7/ A. In each regime, the es- 
sential feature only depends on the weights c? of the ini- 
tial states, and independent of the system-bath coupling 
strength A. In other words, in the framework of the RWA, 
all essential properties for the physical observable are in- 
dependent of the system-bath coupling strength, which 
is not always true with the increasing coupling realized 
in many recent experiments ^45i] . Hence, it is necessary 
to go beyond the RWA. 

A. Exact dynamics and comparisons 

The dynamics of the pairwise concurrence of the two 
qubits in their own baths without the RWA can be ob- 
tained exactly by the hierarchy approach outlined above. 
Initiated from the anti-correlated Bell state ([9]) with 
a = l/\/2, the time evolution of the concurrence from 
weak, intermediate, to strong coupling cases are pre- 
sented in the upper panel of Fig. [TJ The concurrence 
under the RWA has been investigated exactly in Ref. f2^ 
both in the Markovian and non-Markovian regimes. The 
relevant results for the same parameters are also repro- 
duced in the upper panel of Fig. [1] for comparison. 

In the weak system-bath coupling (A = 0.02), the 
rotating-wave term dominates the entanglement evolu- 
tion as indicated in Fig. [TJa), and the results includ- 
ing the CRT show negligible deviation from those in 
RWA. In the non-Markovian regime, the concurrence ex- 
hibits periodic oscillations with amplitude damping. The 
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corresponding time of the zero points is tn = [mr — 
arctan(i?/27)]/i? (n = 1,2,---). Whether it only van- 
ishes at these discrete moments is very subtle, the peri- 
ods of the time for zero entanglement is too short to be 
visible in this case . While in the Markovian regime, the 
concurrence decays exponentially and vanish asymptoti- 
cally. The decay rate is larger than the non-Markovian 
one. So at such a weak coupling, we have not found 
any evident effect from the CRT, so the previous RWA 
description is really valid. 

What happens if we increase the system-bath coupling 
to the intermediate regime, say A = 0.5. It is interesting 
to observe in Fig. [Ub) for A = 0.5 that the dynamic 
behavior is evidently different from that in the RWA. 
From Eq. (11), one can see that if initiated from anti- 
correlated Bell state, in the RWA the ESD never happens, 
which was first found in Ref. [1^ in the same system. But 
without the RWA, in both the Markovian (/ — 5) and 
the non-Markovian (/ — 0.1) regimes, the concurrence 
exhibits ESD obviously. It follows that the RWA can 
not describe the essential feature of the evolution of the 
entanglement beyond the weak coupling regime and is 
therefore broken down under strong coupling. 

Importantly, this physical condition can be realized in 
the recent strong coupling experiments ,45] . In the non- 
Markovian limit (7/A^l), the spectral function of baths 
is reduced to the single mode case by J{uj)fv^S{u! — ujo). 
Hence, the system can be approximately simplified to two 
independent Rabi models. The corresponding effective 

system-bath coupling constant is given hy g = -^/^j so 

g/ojQ — 0.11 for / ~ 0.1, which can be practically realized 
in the recent experiments of the superconducting qubits 
coupled to LC resonators, where the coupling constant 
g/ujQ already reached to ten percentages |45| . 

In the Markovian regime, the entanglement vanishes 
quickly and permanently beyond the weak coupling as 
shown in Fig.[T](b) and[IJc), mainly due to the generation 
of extra phonons with increasing coupling, which disturb 
the pairwise entanglement. No revival of the entangle- 
ment appears in this case. Interestingly, the evolution 
behavior is almost unchanged as the increasing coupling 
strength for A > 0.5, which follows that A = 0.5 is strong 
enough in this case. However, in the non-Markovian 
regime, the entanglement dynamics show essential differ- 
ent behavior with the system-bath coupling. The revival 
of the entanglement after the ESD occurs at A = 0.5, 
but this revival phenomena never happens in the strong 
coupling regime, as shown in Fig. [TJc) for A = 2. We 
propose that the more extra phonons activated by the 
CRT would suppress the feedback from the bath to the 
qubit psl ] even in the non-Markovian regime, resulting in 
the permanent disentanglement. The effective system- 
bath coupling strength for / = 0.1 in Fig. [ijc) can be 
estimated as q/u)Q = 0.45, which is hopefully realized in 
the near future (3,113]. 

Then we investigate the dynamics of the QD for two 
qubits in two independent baths to explore the influence 




FIG. 2: Effect of system bath interaction on dynamics of 
quantum correlation of two qubits of two independent baths 
with initial anti-correlated Bell state 1'1'(-^))- 1^0 = 1, 7 = 
O.IA. 



of the CRT. The QD in the RWA has been discussed by 
Wang et al. |33] and Fanchini et al. \4^. In the lower 
panel of Fig.[TJd), we display the evolution of the QD for 
the same parameters as the upper panel by both the hier- 
archy approach and RWA. In the weak coupling regime, 
it behaves similar to the concurrence in Fig.[lja) for both 
the Markovian and non-Markovian evolutions, where the 
CRT can be ignored. Beyond the weak coupling regime, 
as shown in Fig. [TJe) andlHf), the CRT drives the QD 
to deviate from those in RWA obviously. Specifically for 
the non-Markovian case, the QD decays dramatically and 
shows weak revival signal, which mainly attributes to the 
excitation of the phonons in baths. For the Markovian 
case, the QD decreases faster than that in RWA and de- 
cay to zero asymptotically. 

To see the comprehensive effect of the system-bath in- 
teraction, we plot the dynamics of the quantum corre- 
lation as a function of both t and A in Fig. [2] for initial 
anti-correlated Bell state. Both the entanglement and the 
discord are suppressed monotonically with the increasing 
interaction. The revival phenomena for entanglement is 
not present completely as the coupling strength exceeds 
a critical value, e.g. Ac is about 1.5 for the parameters 
in Fig. [2] (a). The region of ESD becomes wider with 
the increase of the system-bath coupling. But Fig. [2] (b) 
reveals that the QD never vanishes suddenly, consistent 
with the previous observation concluded from an arbi- 
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FIG. 3: Non-Markovian dynamics of quantum correlation 
of two qubits of two independent baths with initial anti- 
correlated Bell states \^{a)) (upper panel) and correlated Bell 
states I^'(q)) (low panel). The insets in low panel are the en- 
larged views, (jjo — 1, X — 2, and 7 = O.IA. 



trary Markovian evolution [5T|. Therefore, such a limi- 
tation is removed in the present observation in the exact 
evolution which definitely includes the non- Markovian 
effect. No sudden death of the QD is universal in this 
sense, in sharp contrast with the entanglement. More 
over, in our exact solution for the full model without the 
RWA, we actually have not found the real zero QD at 
any time, despite arbitrary small. This nonzero nature 
for the QD maybe intrinsic at least for the spin-boson 
model. Previous observation in the RWA that the QD 
can vanish at some discreet times [H, may be an 
artificial result due to the absence of the CRT. 

The essential feature for the dynamics of the quantum 
entanglement and discord under the RWA is not altered 
qualitatively with the system-bath coupling in the same 
non-Markovian or Markovian regime. However, it is not 
that case in the real exact solution without the RWA. For 
the weak system-bath interaction, we have shown above 
that the RWA can basically give the right description 
of the quantum correlation of two separate spin-boson 
models for both Markovian and non-Markovian regimes. 
While in the strong coupling regime, the CRT should be 
necessarily included for correctly depicting the phonons 
generation from the baths. And the quantum correlation 
shows significantly different behavior from that in RWA, 
especially in non-Markovian case. Hence, we focus on 
the dynamics of the quantum correlation in the strong 
coupling and non-Markovian regime in the following with 
various mostly used initial states. 




FIG. 4: Non-Markovian dynamics of quantum correlation 
of two qubits of the two independent baths with extended 
Werner like initial state pAs(O) = r\^}{^\ + \0 = 

|<l>(-i=)) (upper panel) and 1^) = \^{-^)) (lower panel) . The 
insets in the lower panel are the corresponding enlarged views, 
ojo = 1, A = 2, and 7 = O.IA. 



B. Non-Markovian dynamics at strong coupling 

First, we examine the dependence of the dynamics of 
the quantum correlation on the value of a for the ini- 
tial Bell state in Eq. ([9]), which characterize weights of 
the two superposed anti-correlated spin states. Since the 
quantum correlation is similar in the regime a^G(0, 1/2) 
to that in (1/2, 1), the first zone is selected to study the 
correlation evolution in Fig.[3l[a)-[3l^c). In RWA, both the 
concurrence and the QD show damping oscillations with 
stronger amplitude as increases. The other evident 
feature is that they only vanish at some discrete mo- 
ments. When the CRT is included, the concurrence ex- 
hibits BSD after finite time evolution in the whole regime. 
The starting time of ESD is much earlier than first van- 
ishing moment in the RWA. The QD first decreases al- 
most to zero and then revives with amplitude much too 
weaker than that in the RWA, mainly attributed to the 
suppress from the phonons. 

Fig. [2Id)-|3Uf ) presents the results for the another ini- 
tial Bell state in Eq. (fTO|). It is known that the discord 
dynamics for Q:^e(0, 1) is similar to the concurrence for 
Q!'^g(1/2, 1). Hence, we also select these initial states for 
comparison as in Ref. (43l |. The exact dynamics shows 
that the ESD occurs very quickly and irrelevant with 
the value of a. The QD first decreases almost to zero, 
and then revives weakly. The revival amplitude becomes 
stronger with the value of a. For the same a, the revival 
amplitude is much too weaker than that in the RWA, due 
to the CRT effect, as demonstrated in the inset. 
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Finally, we assume the extended Werner like state 
{pab{0) — i^lOm + "'t^I) the initial state and see 
what happen to the evolution of these quantum corre- 
lations. The main results are list in Fig. S] As one 
classical kind of the mixed states, it is important in 
quantum information processing |52h54|. As the initial 
state with |^) = \^{^)): the concurrence in RWA shows 
complete death in re(0, 1/2), whereas the discord in the 
regime r€(0, 1) behaves similar with the concurrence in 
r€(l/2,l) [33]. Hence, we study the quantum correla- 
tion in the regime rS(l/2, 1) to explore the effect of the 
CRT in Fig. |ll[a)-|3^c). The entanglement vanishes sud- 
denly and permanently for all values of r, RWA ones 
show damping oscillations for large r (e.g. r = 0.7, 0.9). 
The QD is also strongly suppressed first, and then revives 
weakly. The revival amplitude increase slightly with r. 
It is again shown that the QD is not so fragile as the 
entanglement. If |5'(--^)) is replaced by |^) in the initial 
extended Werner like state, no essential different behav- 
ior can be observed except that the revival amplitude 
becomes much too weaker, as demonstrated in Fig. IDJd)- 

For the three kinds of initial states studied above, it is 
generally find that the larger initial quantum correlations 
is, the stronger they evolves at the same moments. 



IV. TWO QUBITS COUPLED TO ONE 
COMMON BOSONIC BATH 

The Hamiltonian of two qubits coupled to the common 
bath H = Hg + Hh + Hsb reads 



^« = E 



■i=A,B 



k 

Hsb = Vj2{gkh+glbl), 



(13) 



where the notations are the same as those in Eq. (2). If 
the spectrum of the common bath is also Lorentz distri- 

■, and the correlation func- 



2tt {uj-uJoV+l' 



bution J{uj) 

tion of the bath should be also C{t) = ^e-(^+*'^'')*. By 
choosing the proper auxiliary influence functional in Ap- 
pendix, the hierarchy equation is given by (|A17p (4l| 



dpn{t) 
dt 



= -{iH^ + n ■ v)pn{t) ~ i y'pfi+eAt) 

k=l,2 

fc=l,2 

(14) 



where n = (ni,n2), the unit vector ei = (1,0) and 62 
(0, 1). The reduced density matrix for qubits is pAB{i) 



P(o,o) (Oj by which we can derive the quantum correlation 
numerically. 

In RWA, the system-bath interaction becomes Hf^^ — 



+ ^+) Efe 9kbk + {o^ + o^) Ek albl If we repre- 



sent the Hamiltonian in the basis {|0) — |00), | -f ) = 
i(|10) + 101)), I - ) = ^(|10) - 101)), |2) = 111)}, I - ) 
is found to be isolate from others. Hence we can rewrite 
the Hamiltonian in three level formlSSi 



Ho = 2tJo|2)(2|+wo|+)( + |-c^o, 
H,nt = V2^5fcafc(|+)(0| + |2>( + | 



(15) 



H.c. 



For the single particle excitation or double excitation 
case, the motion equation of the system density matrix 
can be derived by the called pscudomode mater equation 
in the interacting picture [43 . ,55 - .59j 



dpijt) 
dt 



= -^[V,pi{t)] 

--f{a'^api{t) + pi{t)a)a - 2api{t)a^)ilQ) 



where (a) is the creator (annihilator) of the bosonic 
pseudomode, 

y= V^(a|+)(0|+at|0)( + |+a|2)( + | + at|+)(2|). 



and pi{t) — 



p{t)e 



-iHnt 



By integrating the bosonic 



part of the density matrix, the system reduced density is 
given by pAB{t) = Tra{/o(i)}. For single excitation, an al- 
ternatively exact solution has also been done in Ref. [slj . 
If the initial state is chosen by X form in Eq. ([6]), the 
concurrence is derived by [59| 

C{t) = 2max{0, IP23I - \/piimJpi4| - \/P22P33}- 

(17) 

For the common bath, the entanglement evolution ini- 
tiated from the anti-correlated Bell state without the 
RWA has been studied by Ma et al. 41]. On the other 
hand, in the previous studies of the dynamics of QD un- 
der the RWA, the correlated Bell state is usually selected 
as the initial state [H, [H, Isl, [S^. Therefore in the 
present exact study without the RWA, we will focus on 
the initial correlated Bell state. 

First, we exhibit the dynamics of the pair wise concur- 
rence of the two qubits coupled to a common bath with 
and without the RWA initiated from the Bell state in 
Eq. © with a = 1/ -\/3 in the upper panel of Fig. [5] for 
various coupling strength in both Markovian and non- 
Markovian regimes. Similarly, for the weak coupling 
(A = 0.02), the CRTs take little effect for both dynam- 
ics in different evolution regimes, showing the same re- 
sults as in Ref. [i^. In the intermediate coupling regime 
(A = 0.5), the ESD occurs more quickly than that in 
the RWA in the Markovian regime, and the evolution 
shows slight deviations from the RWA ones with damp- 
ing oscillations in the non-Markovian regime. When the 
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FIG. 5: Dynamics of quantum correlation of two qubits cou- 
pled to the common bath with initial correlated Bell state: 
|*(;i=)>. .^0 = 1, 7 = /A. 



FIG. 6: Dynamics of quantum correlation of two qubits cou- 
pled to the common bath with initial correlated Bell state 
|vl'(-l=)) (a) concurrence and (b) QD. ixjq = 1, 7 = O.IA. 



strong coupling regime is reached, the CRT takes remark- 
able effects in both evolution regimes, where concurrence 
vanishes completely and permanently without the RWA. 
Hence, the CRT also suppress the entanglement dramat- 
ically as the system-bath coupling for the common bath, 
similar to those for the independent baths. 

Then we turn to the relevant QD, as shown in the 
lower panel of Fig. [5] It is interesting to note that in 
the non-Markovian regime, the CRTs take very limited 
effect with the increasing coupling. Even at the strong 
coupling, only slightly deviations are observed. More in- 
terestingly, the QD takes finite value at any time, and 
take a moderate finite value asymptotically in the long 
time limit, for arbitrary coupling strength. In the Marko- 
vian regime, the QD tends to a remarkable value in the 
intermediate and strong coupling regime. Actually, even 
at weak coupling like A = 0.05, the asymptotical value of 
QD takes a considerable non-zero value. 

From the numerical iteration of the hierarchy equation 
in Eq. (ITi| with Bell correlated initial state under Markov 
regime (e.g. / = 10 in Fig. [5]), we find that the form of 
the two qubits reduced density matrix under the same 
basis as in Eq. (|6|) always to be X type as 



Pit) 



a{t) w{t)\ 

b{t) b{t) 

b{t) b{t) 

Kw*{t) d{t), 



(18) 



In the weak system-bath coupling (say A = 0.02), a{t) 
and b{t) decrease gradually to after the finite time 
evolution, whereas d{t) rises to reach 1. Hence both 
concurrence and the discord decay to exhibit mono- 
tonically as indicated in Fig. [SJa) and[5Kd), similar to 



that in the RWA. When the coupling strength becomes 
stronger, e.g. A — 0.5, 2.0, the phonons are excited con- 
siderably due to the strong coupling [i^. According 

to the CRT 'Ek9*k^li^+ + 0-+)' the qubits are also ex- 
cited accompanying the extra excitations of the phonons, 
which stabilize a{t) and b{t) in long time evolution. 
It is interesting to observe that they tend to 1/3 and 
1/6 in the strong coupling limit at t— >oo. The corre- 
sponding reduced density matrix is then described by 
p{t^oo) = i(|00)(00|-h|ll)(ll|-h|-h)(-h|). which is not 
altered by different a^. Under this separable mixed state, 
the disentanglement occurs naturally from Eq. (jl7p . and 
the QD is stabilized at 0.23. 

To see the overall effect of the system-bath interac- 
tion on the dynamics of the quantum correlation beyond 
the weak coupling regime, we draw the corresponding 
3D plots for A > 0.1 in Fig. [5] in the non-Markovian 
regime. For the entanglement, one can find that there 
are two regimes. One is damping oscillation at A < 1. 
The other is the regime where the entanglement revives 
following the ESD, and the revival will be suppressed and 
finally disappear with the further increase of the cou- 
pling regime. Very interestingly, after weak oscillations, 
the QD is generally robust in the long time evolution 
in arbitrary coupling regime, which is very useful as the 
quantum information resource. 



V. CONCLUSION 

In summary, without performing any approximation, 
we exactly investigate quantum correlations in two typ- 
ical kinds of spin-boson models by applying the hierar- 
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chy approach. The corresponding RWA results are also 
reproduced for comparison. The Markovian and non- 
Markovian effects of two kinds baths in strong system- 
bath interaction on the dynamics of the quantum cor- 
relation are revealed in detail, which demonstrates that 
the CRT should be necessarily included beyond the weak 
coupling regime to capture the essential features of quan- 
tum correlations correctly. 

For two independent baths, the CRT suppresses both 
the entanglement and the QD dramatically, due to the ac- 
tivations of phonons in baths as the system-bath coupling 
increases. If initiated from the anti-correlated Bell state, 
The ESD is driven by extra phonons activated with the 
increasing system-bath coupling, which is absent in the 
RWA study. It follows that previous picture under the 
RWA is essentially modified with the increasing coupling. 
The QD is found to be always higher than zero, despite 
extremely small, also different from that in the RWA. The 
dynamics with other initial conditions is also discussed at 
the strong coupling and non- Markovian regime. The gen- 
eral trend is that the larger initial quantum correlations, 
the stronger they evolve. 

For the model of two qubits coupled to the common 
bath, the suppression of the entanglement by the CRT 
is also observed. It is interesting to find that the non- 
classical correlation reserved by the QD can be enhanced 
as the coupling strength, in contrast with the entan- 
glement. It follows that while the entanglement be- 
comes fragile with the increasing coupling to the sur- 
rounding environment, the QD is however robust in the 
long time evolution. In recently experiments, the strong 
system-bath coupling has been reached in circuit QED 
systems [i^, where g/a;o~0.1 for multi modes and the 
RWA description is broken down. The present study 
for the evolution of the quantum correlations beyond 
the RWA may motivate the corresponding experimental 
studies based on these strong coupling systems. 
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Appendix A: basic concept and application of the 
Hierarchy equation 

The Hamiltonian of the system-bath coupling system 
is modeled a.s H = Hg + Hi, + Hgt ■ Then the time de- 



pendent density matrix is given by p{t) 



-iHt 



p(0)e 



e ^ lo Under the specific representation a)}, 

{xt,at\pi{t)\xt,at) 
= {xt,at\Ui{t)p{0)uj{t)\x't,a't) 

= J dxodxQ J daodaQ{xt,at\Ui{t)\xo,ao) 

X {xo,ao\p{0)\xQ,af^}{xQao\Uj{t)\xt,at). (Al) 

It is known that pl{t) — Tr(,{p/(t)}. Hence by integrating 
the bath degree of freedom, the reduced system density 
is 

{atWsit)\a't) 

dxt{xt,at\pi{t)\xt,af.) 

= J dxtdxodxQ J daodaa{xt,at\Ui{t)\xo,ao) 
X {xQ,aQ\uj{t)\xt,af.){xQ,aQ\p{0)\xo,aQ). (A2) 
Considering the propagator in functional form 

{xt,at\Ui{t)\xo,ao)= H '^i^i^)] 2?[a(T)]e^^[-M'"M] 



and the initial factorized state p(0) = Ps{0)^Pb [pb = 
exp(—f3Hb)/Z is the equilibrium state), the reduced den- 
sity matrix is reexpressed by 



{at\psit)\a[) 
= / dandaQ{an\ps{t)\aQ) 



(A3) 



at 



at 



In the interacting picture, the corresponding density 
matrix reads pi(t) = Ui{t)p{Q)U]{t), with Ui{t) = 



X / 'D[a{T)] I ^ V[a (t)] / dxtdxodxQ{xo\pB\xQ) 
X / ' V[x{t)] f '2?*[a;'(T)]e*^[^(^)'"(^)l-'^[^'(^)^"'(^)l. 



Then the functional field is introduced by 

T[a{T),a' (t)] 

dxtdxodxQ{xo\pb\xQ) 

X / ' V[x{t)] I 'l?*[x'(T)]e*^[='(^)'"(^)l-''^[^'(^)'"'(^)l 

J Xq J x[, 

= Tn{pbUl[a]Ub[a]}, (A4) 

where {xt\Ub[a{T)]\xQ) = /^""J X'[a;(T)]e*^[^(^)'"(^)l. The 
functional field generally satisfies J^[a, a]<l, where 
J^[a, a ] — 1 if there is no system-bath interaction. 
Here, we include the harmonic bath with the Hamil- 

tonian as Hb — J2k 2^ + '"^ xl- Without loss of 
generality, we firstly consider only one mode case as 
Hb " + ^^^i^, and the system-bath coupling term 
is Hsb — Qx. In the interaction picture, Hsb = Q{t)x{t), 
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and Ullair)] = g-'/o Through the tedious 
by standard Feynman- Vernon influence functional calcu- 
lation, we derive the influence functional as J-[a,a ] = 
g-$[Q,a ]^ where the influence phase is 

= [ dt I ds{Q[a{t)]-Q[a{t )]} (A5) 
Jo Jo 

X {Cit - s)Q[a{s)] - C*{t - s)Q[a\s)]}, 
with corresponding bath correlation function 

C(r) = ^[coth(/317/2)cosr2T-isinl7r]. (A6) 
mil 

Here we should note that the variable Q[a] is a real 
number. Then we generalize the results to the multi- 
mode case, where the system-bath interaction now is 
Hsb = J2a QaFa, and Fa = J2k '^k^k- The influence func- 
tional becomes 

^[a{T),a{T)] = exp(-^ / dr{Q,[a(r)] - Q^a' (r)]} 

a Jo 

x{Q4a{T)]^Q:[a{T)]}), (A7) 

where 

QaHt)] = V / drCaa'it - T)Qa'[aiT)]. (A8) 

a' 

The correlation function in independent baths is 

Caa'ir) - {Fa{T)Fa,m (A9) 

= E|S&['^°th(^)cos(l^^)-^sin(17^)]. 



The proper way is employing the hierarchy equation by 
introducing the auxiliary influence functional 



Ps{n,t) = / (iaodao(ao|/Os(0)|ao) 



(A12) 



V[a{T)] / V*[a (T)]J-„[a,a']. 



While such auxiliary choice is closely dependent on the 
specific form of the bath correlation function. The hier- 
archy equation used in the present paper is described in 
the following. 

For two independent baths with Lorentz type spectral 
function Jq(cl>) = -^-i — ^ili , — i, 
simply is 



the correlation function 



\ala 



exp[-(7, + iwg)^]. (A13) 



By using the spectral function Ja(w) = | X]/c man° ^{'^ ^ 

k k 

ri^), we finally have 

(iwJa(w) [coth(-|-) cos(a;r) - isin(a;r)] 

In the schrodinger picture, the motion equation of the 
system density is 



according to the following formula 

we introduce the auxiliary influence functional of each 
subsystem as 

Fai[a,a] = (Q4a(t)])'=^(Q:[a (t)])'==J-[«,a'],(A14) 



where Qa — cr^, (a — A, B) and index k — (fci, fc2). Then 
the whole influence functional becomes 

Ffijri[a,a ] =FA,n[a,a ]FB,rh[a,a ]. (A15) 

Finally, the evolution of the density matrix reads 

d 

-Q-^Pn,m{t) = -{iHg +n- UA+m- VB)pn,fa{t) (Al6) 



k=l,2 



n,[Q'X + (,-lf+'Q'^]pn-e,Mt) 



dpsit) 



iHst JHst 
dt ' "L"0,rovvj ' - dt 



-i[Hs,ps{t)]+e- 



(AlO) 



k=l,2 

fc=l,2 

Ab7b 



From Eq. (|A3|) and (|A4)) . we have 
dp's^t) _ 



A;=1.2 



dt 



daodQ;o(ao|ps(0)|ao) (All) 



"0 



V[a{T)] , V*[a{T)]- 



dt 



The key step here is to connect time differentiation of 
the density matrix with that of the influence functional. 



where unit vector ei ~ (I7O) and €2 — (0,1) Va ~ 
with v'^ = 7a + and vt = 7a — «i^o- 
A^B = AB ~ BA and A°B = AB + BA. 

For common bosonic bath with the similar Lorentz dis- 
tribution, if we choose the auxiliary influence functional 

^n[a, a] = (Q* [a (t)])"^ (-j)"^ (g[a(t)])"^ J-[a, a], 
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the hierarchy equation can be expressed by 

fe=l,2 



4 

fe=l,2 



(A17) 
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